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ABSTRACT 

Recently many studies have considered the possibihty of a Lorentz Invariance Violation (LIV), and ex- 
plored its consequences in a wide range of experiments. If this is true, a LIV could explains some mysteries 
in Cosmology. In this paper specifically, we will analyze the effects on The Primordial Baryogenesis be- 
cause it is one of the more important and mysterious phenomena of the Big-Bang, that happened at very 
high energies, so we have a real chance to obtain an important effect. We will see that this effect could 
exist, depending directly on the temperature, that is very high at this time in the history of the Universe. 
So, it is possible to use this result as a test for a LIV and explore the possibility that the boson that 
started the baryogenesis explains, in part, the dark matter. We will obtain estimates about the beginning 
time of the baryogenesis and the boson mass too, that come directly from the LIV. 

INTRODUCTION 

In 1928 Paul Dirac, in his famous equation, predicted the existence of anti-particles. Now this has been 
experimentally checked: to each particle one anti-particle is associated, where both have the same mass, 
half life and opposite charge. It is natural to think that they must exist in the same number, but the 
reality clearly shows the opposite; it shows an asymmetry baryon-antibaryon. The production of both 
and this asymmetry is named Primordial Baryogenesis. 

What Baryogenesis's theories have in common are the _B, C and CP violation, maintaining CPT 
symmetry, and Departure from Thermal EquilibriurrQ. One of the favorite theories to explain the Baryo- 
genesis consists of introducing a new particle named X Boson that is of the same form as Z and W 
bosons, and should be related to a new force that is necessary for the Grand Unification Theory [2]. 

We will study this phenomenon using statistical mechanics and cosmology, but we will introduce a 
Lorentz Invariance Violation {LIV) too [3]- [11], that is translated in a small modification of the usual 
energy dispersion relation. 



In the Section 1 we will show the usual procedure to represent the Baryogenesis in the most simple 
model, explaining in detail the conditions previously mentioned. In the Section 2 we will show the 
effects when the LIV is included through Threshold Energy conditions and we will incorporate them to 
the development of the previous section. Finally we will expose the results and some predictions obtained 
by the LIV. An extra effect by the LIV is a possible breaking of CPT symmetry, but we will forget it 

^These conditions were proposed for first time by Sakharov [l]. 
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because these effects turned out to be very little. 

The LIV that, apparently, is the most important at high energies, and that we will use here, has the 
form [6] [2] [9] [lO] [n]: 



E^^vl,y + m,'c'' (1) 

where v^ax = c(l — a) is the maximum particle velocity, a is very little such that it will be marked 
to high p. Some estimates of this factor a, given by theoretical analysis of the physic of each particle and 
experimental limit of some parameters, say that is approximately 10^^^ or 10^^'^. This value will always 
be little but it can be different for distinct particles [7] [TT] . 



1 BARYOGENESIS. 

For Baryogenesis to be possible, we need a Boson X that produces baryons and anti-baryons when 
decaying. Three conditions have to be also carried out that will allow the present baryonic asymmetry 

a) Baryon number violation. This mean that, while the Universe expansion is happening, the factor: 



where nb, ni and s are the baryon and anti-baryon numbers and the entropy per comobile volume, 
must change. If B violation had not existed, the asymmetry would be only given due to the initial 
condition because B = cte. So, if we consider a, X = X model, there must been a reaction of the kind: 



X < — >b + b 

X< — >b + b (3) 

b) C and CP Violation, with CPT symmetry. Of this form baryons and anti-baryons have the same 
physics (same mass, Lorentz violation, etc) and that the boson decay rate can be different (in a small 
factor) to the inverse boson decay rate. If these violations did not exist, baryons and anti-baryons will be 
created in the same number, without preference for the matter above the anti-matter. We can interpret 
the CP violation, in term of the amplitude of probability, as: 



I M{X ^ 6 + 6) p - I M{b + b^X) p= ei 

I M{b + b^X)\^ -\ M{X -^b + b) p== 62 (4) 

where 1 ^ > allow a little inclination in favor of the matter. As we also want to keep CPT, we 
must have: 



M{X ^b + b) p = | M{b + 6 ^ X) p 

\M{X ^b + b)\^=\M{b + b^X)\^ (5) 



So, we can deduce the most simple model: 
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I M{X -^b + b) p=| M{b + b^X) |2= 1(1 + e) I Mo P 

I M{X^b + b) Mib + b^X) |2= 1(1 -e) I Mo |' (6) 
Where | Mq P is constant. 

c) Departure from Thermal Equilibrium. If equilibrium exists, we permanently have that Ub = n^, 
and the matter and anti-matter amount would always be the same. A Departure from Equihbrium 
process is described by the Bohzmann equation : 

L[/] - C[/] (7) 

Where L is the Liouville operator, C is the Cohision operator and / is the distribution function. The 
Liouville operator corresponds to: 

and the Collision operator must consider all possible combinations of the process. Evaluating: 

Using the FRW metric and considering that — [E.VmaxI^ and a;^ — [vmaxt^x], is reduced to: 

E df R _ , ,2^/ 



Applying / it is obtained: 



VmaM'^ = C[/] (10) 



9 f 9f ,3 rju^^-Vraax fp'^^f 3 gVmax f^,^^<Pp 



d'p^Hit)^ C[f]^ (11) 



{27rh)^ J dt " {2TThf J EdE (27rft)3 J E 

Where H{t) = is the Hubble constant. Integrating by parts the second term on the left and using 
the definition of particles number density n(t) = ^2nh)'-^ I f'^'^P^ obtain: 

^,3.,«„,«.|^/Cl;4 (12, 

Where the term of the right represents the departure from equilibrium. If it is declared null, it is 
obtained n{t) =^ nQR^^{t) that is the case of a evolution in equilibrium. 
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2 THRESHOLD ENERGY AND COLLISION FACTOR. 

As in this paper we want to see Lorentz Violation effect, we can have energy levels where the reaction is 
not produced. To obtain this energy level, we will use the Threshold Energy [5] [10] [T^], that consists on 
the following: 

The total energy of a system of many particles is: 



E = J2E.d\p.\)+^,i^P>-Y.pij (13) 

where are Lagrange multipliers that impose the conservation of momentum, i level to every particle 
and j represent the vectorial components. Deriving with respect to p], to minimize E, we obtain: 



S^-HO (14) 
dpi 



So, E is minimum when the velocity is the same for all particles. We will name this energy Emin = 
Ei where all particles i have the same velocity. Now, if we have two particles that collide to produce 
others, we can write the energy as: 

E = E, {p- \p2\n) + E2 {\P2\) + x{n^- l) (15) 

where fi is normalized vector of p2 , P = Pi + P2 is the total momentum and x is other Lagrange 
multipliers that impose the normalization of n. Maximizing E with respect to n, we obtain: 



2x 

Since |7t-| = 1, we have that x = ±-^^^-j^, so: 

n = ±^ (17) 

Evaluating in E, we can see that the maximization occurs when n = — — , so when a frontal collision 
happen. This case is represented for E^ax = Ea+ Ef, where a and b are particles that collide face to face. 

Finally, the threshold condition is given by: 

Ejnax ^ Ejriin (1^) 

Since a and b go from one to the other and the particles i go in the same direction, we have: 
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2.1 Reactions Allowed Zones 
2.1.1 X — > 61+62 

hi can be a baryon or anti-baryon. Using (jl9p : 

Ex > -Ebi + Ei,2 = 2£'f, 
Px = Pbi + Pfcs = 2pfc (20) 

as they are the same kind particle and have the same velocity, we have that pb^ — p&j = pt and 
= = Eh- As a dispersion relation is carried out of the form E^ — WmaxP^ + m^c^, we have: 

- v'x)p\ < {m'x - ^mDc" (21) 

If Vb > Vx- 



Vx<. rV^f\ c- V.<s\ '2~^7/ (22) 

V K - ^x) V 4(^6 - ^x) 

giving an important superior restriction to momentum. If we use mx 3> Wb and v\ — v\ = [vf, + 
vx)ivb — Vx) — 2c^da con da — ax — ctfo, it is reduced to: 

mxc mxc 

- - ^^^^ 

V2aa 2v2oa 

If Vb < Vx, we do not have a bound, since p^, > 0. 
2.1.2 61 + 62 — >X 

Doing the same development using we obtain: 

Ebi + Eb2 > Ex 
Pbi - Pb2 = VX (24) 

AvIpI^ - AvlpxPb, + p\ [vl -v\)- rn'xc^ > (25) 

Where the approximation mx 3> to& was used. The solution, in pb^ , is a parabola with a minimum. 
So the zeros will give us the bounds of the reaction. If we define f{pbi ) = ap^i + ^Pbi + we will see that 
a = Av^, b = —Av^px and c = Px{vb — Vx) — ra\c^ where the zeros are given by: 

Pb^ ,0 = A2 = 6^ - 4ac 26 

la 

If < 0, zeros do not exist and f{pb) > is always carried out without bound. On the other hand, 

2 7712 
b^x 



if A^ > 0, a zone exists where the reaction is prohibited. Evaluating, it is seen that A^ = IGv^Ej^ > 



that is, a bound exists. The zeros are: 



Vbpx ± Ex , . 

Pbifi^ 7y 27 



So, the bounds for this reaction are: 
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P"^- 2^, P'^- 2^, ^2^) 

As and are related to pb^ — Pb^ = Px , that represent the two particles when directly collide 
face to face, so sign{pb^) = signljph^)- Comparing: 

For the first bounds: 

VbPx + y'v\p'x + "ni^ijC^ -VbPx + Vv'xPx + "^x^^ .^q^, 

Pbi > 7. Pb2 > 7, (^y) 

2vb 2vb 

If Px > 0, clearly always pb^ > 0, but pi,^ > will only be if {vf — v\)p\ < m^c^. If px < 0, pb2 > 
and pbi > only if {vf — v\)p\ < m\c^. 

For the second bound: 



VbPx - V^jPx + m"]^ _ -VbPx - \J'v\p\ + 



2 c4 



Pbi < 7 P&2 = -Pfci - < ^ (30 

2vb 2vb 

In this case, if px > 0, always pf,^ — Oi but pb^ < will only be if (u^ — Wx)^^: 5; mj^c'^. The same 
happens if < 0. So, the condition to signijpb^) — sign(jib2) is that: 

{vl - vDpI < va\c^ (31) 

That corresponds to the decay condition that in X — > hi +62- Analyzing, we can simplify ((28|) . 
considering \px\ = \Pb^,A - \Pb^iX in: 

I I . Vb\px\ +Ex I I -Vb\px\ + Ex , . 

\Pb^^^ I > 7, I > 7, (32) 

2vb 2vb 

but pb^ can be in any of both regions. So, the bound can be simplified to: 

^ -vb\px \ +Ex , . 

- 2^b ^^^^ 

That, together to ([3T|) . represent the allowed zone for the reaction. Now that we have the thresholds, 
we can do our calculations. 



2.2 Collision Factor and Departure from Equilibrium Condition 

The Collision Factor is: 



{2nhy 



■ Jc[f]^ = - j {2TTh)H\px -Pb.-Pb,) TxMM dUidU2dIlx (34) 



TxMM = /x (I Af (X ^ 61 + 62) I' + I il/(X ^ 61 + 62) H 

-AJfc. I M{bi + 62 -> X) |2 -/jj^^ I M(6i +h^X) |2 

Where fx, fbi and fi. are Boson, baryons and anti-baryons distribution functions respectively, and 

'^^^ ~ (2^1)^ '^ET' ^^'^ ~ (^^ffi ■ "^^^ amplitudes are given by the model previously mentioned 
but are cancelled out of zone of calculated thresholds. This expression follows from assuming that the 
distribution of boson, baryons and anti-baryons approach a Boltzmann distribution of the kind: 
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fx =e 
hi = e' 

4 ^e-^ (35) 

Where ^x, M =''^6 the boson and baryon chemical potential. We use this to simplify the expression 
of the collision factor. This is acceptable in the high temperature approximation, that we are using 
throughout this work. 

We must mention that the boson is not decoupled, since if this is not like that, the effect on the baryonic 
asymmetry would be little for being practically Non-Relativistic (T < mx)- This is that bosons, baryons 
and anti-baryons are still in chemical equilibrium with the thermal bath. So we have that = —^J■b■ = ^J■■ 

Analyzing the product fbifb2-, we can see that: 



2m 

k /3 k dT 



e e 

req 

Ix 



/fe^ (36) 



Where is the chemical equilibrium boson distribution (/ix — 0). The same way for the product 

fbj-b2=fx'^~'^ (37) 
With these relations and using our probability of amplitudes, we obtain: 



(27rft) 



/ C[/]^ = -J{27rhf I Mo f 5\px - Pb, -PbJ x 



[fx-fx']dTlidIl2dnx (38) 



Where we have considered that fj, ^ A^T and we keep up to first order in and e. As the term fx 
comes from the boson decay and f^ of the inverse decay, they have different integration ranges. For this, 
we must separate them in the integral, so: 



Where I^'^ ^ and /j,^ '^^^ ^ contain the distribution with and without equilibrium respectively and 
Pi = [Ei\ Vmax,iPi]- We wiU rcsolvc them in general form and then we will distinguish them. For this, we 
use the relation [131: 



vl<Ppb2 
2Eb2 



d*Pb2S {pI - mlc^) Q{Eb2) (40) 
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Then we integrate in d'^Pb^ using S'^{px — Pbi ~ Pb-i)^ that is simply replacing p^^ = px — Pb^, having: 



Now we use the second delta to integrate in cpi. For this, we have the identity: 



(41) 



With: 



-f" (</>!) = Px+Pbi - '^PxPb, - ml^c^ 

= m\c^ - 2ExEb^ + 2vbVx\pbA\px\ {cos{0x) cos(6ii) + sin(6'js:) sin(6'i) cos(0i)) 
= -2vbVx\pbA\Px\sin{9x)sm{ei)sm{(t)i) (43) 

We can see that two values of exist; The first between and n, where sin(0i^j) > 0, and other 
between tt and 27r, where sm{(j)i^i) < and equal in module to the previous one. So, we evaluate in 
(<^i i)l with one of them and multiply by 2, obtaining: 



oetween tt and Zn, where sm(0i_ij < U and equal m m( 
P'{<l^i,i)\ with one of them and multiply by 2, obtaining: 

W0i)) = jF7^^^i<Pi-M(^{F''i^i,i)) 



IwM,x = 2 j ^^^^^^^ Ex\F^{i!'^)\'^'''^^ ^^'^ ''l^^-^l' ''"(^^^ sin{Ox)dpb,derdpxdexdcl>x (44) 
Where 6 (-F'^((/>i,i)) appears to assure that cos^(^i) < 1. Besides, the value of carries out: 

2{ExEb, - Vbvx\pbi\\px\ (cos(6'jf)cos(6'i) +sin(6'x)sin(6'i)cos((/)i,i))) = mjfc"* (45) 
Therefore: 

|F'(</.i,i)| = Vacos2(^i) + 6cos(^i)+c (46) 

With: 

a = -{2vbVx\pbi\\px\f 

b = -4vbVx\Pbi\\px\cos{ex){-2ExEb,+mj^c^) 

c = {2vbVx\Pb^\\px\sm{ex)f -m%c^ -AE%El+4ExEb,m%c' (47) 
So, the integral in 6i is reduced to: 

r ^i-£±t±Sd. (48) 

j„i V ax^ + &a; + c 

With X = cos(0i). The parabola ox^ + + c has a maximum (a < 0) and the zeros are within —1 and 
1. This mean that the integral interval can be extended to [—00,00] without affecting anything thanks to 
the Heaviside's 9. Using the relation: 
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) -dx = ^=6(6^ - 4ac) 

y/ ax'^ + bx + c \/—a 



(49) 



we obtain: 



ja 



v^vx 



ExPx 



dEb^ d^px 



(50) 



Where we have used that Ei,-^ = v^pi,^ . Analyzing the second Heaviside's, we can see that its argument 
is positive if: 



Ex - vxpx ^ ^ Ex + vxPx 

o — 01 — o 



(51) 



Now we must distinguish the following processes. 



2.2.1 X Decay 

(a = Neq) From ([23|) . we have the bound: 

mxcvb 



Eb, < 



V px< 



mxc 



Eb^ < oo V px < oo 



If: Vb > vx 



If: Vb < Vx 



(52) 



(53) 



But we also must carry out the limits imposed by the Heaviside in the integral. With a bit of analysis, 
we can see that the bound given by the Threshold Energy is always greater than the inferior Heaviside 
limit. Therefore, the integration limits are: 



Ex- 



vxPx <^ p <^ Ex+vxPx 
2 ^ -C/fti :i 2 



(54) 



If: {vb <vx)o {vb > Vx) and (^px < 



2^2da , 



Ex—vxPx <- p <- mxcVb 



"1 - 2^/20^ 



(55) 



If: 



Where we have used that da ^ 1. Therefore, evaluating in ([5D|) when Vb < vx, we obtain: 

Ibublxi^b <vx)^^J "^d'Px With: < < oo (56) 
These is almost no difference to the case without Lorentz Violation. But, if Vb > vx we have: 



h,M,xi^b > vx) ^ -J- 

"b " 



v^vx 

f fx ,3 , f fx ,3 i f fx 

' / -E;-d Px + yvb / — d px - t; / — 

J A Ex JbExPx 2JbPx 



Vx f fx 



I —d^px + -TT / T^d-'PX 

B Px 2 Js Ex 



(57) 
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with A (px ^ y) and B ^ {y < px < 2y) where y 



We can see that many extra factors 



2V2da ■ 

due to the Lorentz violation and his prohibited energy zones appear. If we call C —> {px > 2?/) and 
consider that in the B and C region have that Ex = vxPx, we obtain: 



lHZ,xi^b >vx) = -3 



fx ,3 , 47rywf, 
—d PX + — 3— 



Where the integration zone in the first integral extend to all momenta. 



47r f 

fxdEx 3- / ExfxdEx 

B '"x Jb+c 



(58) 



2.2.2 Inverse X Decay 

(a — eq) In this case, the bound is given by (|^ and so: 



Eb, > 



-VbPx + Ex 



with: 



(59) 



. mxc 

Px < rrr^ ^b > VX 



V2da 

px < 00 If: Vb < vx 



So, the limits, considering the Heaviside of the integral, are: 

Ex - vxPx ^ ^ Ex + vxPx „ 

< Eb^ < If: Vb > vx 



Ex - VbPx . p , Ex + VxPx „ 

< Eb, < If: Vb < vx 



Therefore, if Vb > fx, we have: 



jeq 

^bi,b2,X 



{vb > Vx) = -3 





r il 




A+B Ex 


TT 




^^6 . 


1 Ex 



d^px 



47r 



"X Jc 



^x 



(60) 

(61) 
(62) 



(63) 



Where we have used that in the region C the boson is ultra relativistic. And if Vb < vx- 

I!lb2,xivb < vx) 



T^{Vb + Vx) f fx ,3 

-^d-'px 



2VbVx J Ex 

req 



nuA-b <-x) = ^J ^d'px 



(64) 



Where we use that ^^^^^ < 1. 

vx 
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2.3 Differential Equation Solution and Analysis 

The distribution functions inside and outside of equilibrium are related by fx = /^^ e'=B^. Now we will 
evaluate in (l39t. 



a) If vt, < vx, we use ([55)) and (IM)) : 



' 3H(t)n(t) = I Mo 



5f 



4u4(27rfi)5 



P (e^ - l) I gdV 



(65) 



b) If Vh > Vx, we use ([55)1 and 



Si 



Ex 



_I^X_ 4:71 

v 



X \ JB JB+C 
eq 



Exfx'dEx 



^ f fx'ExdEx 



(66) 



As the unique part that depends on the time in is the temperature T, if we derive it, we have 
= ~^Ex fx 7 with /3 — So, deriving the differential equation and remembering that — 

I fxd^Px, with nx = e'^'^^n^^, we obtain: 



a) If Vb < Vx- 



nit) + 3 H{t)n{t) + H{t)n{t) 



= M{t) [n''^{t)-nx{t)]+tix 



dp 



dt e/^Mx „ 1 



[n{t) + 3H{t)n{t)] (67) 



b) If Vb > Vx- 



n{t) + 3 H{t)n{t) + H{t)n{t) 



With M{t) = I Mo I 



1^ f and 



M(t) 



n^{t)^nx{t) + 



ingx dJ 



i;|(27rft)3 df3 
+ ^^^;]3^r^['^W+3if(0n(i)] 



9^ e^^^^ 



(68) 



, / f^^^dEx - f Exfx'dEx) + I Exfx'dEx 
Jb Jb+c J Jc 



J = e^''^ I yvb 

It is the factor that represents the Lorentz violation effect. Integrating: 



(69) 



J ^ ^e-^y'"' [e-'^y- {2f3yvb + 1) - e'^'^^ {Pyvbe-^y"" + l)] 



(70) 



By the high temperature, we know that Pyvb ^ 1, moreover if the reactions that produce the Baryo- 
jenesis are sufficiently fast, we have that P/j-x <C 1. So: 
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As If = iiHit), so M{t) = ^g'l^f^, I Mo P > 0. This means that: 

F {fix, fix, nx, fix) « -H{t)T^ cx -T^ (72) 

where F is the usual differential equation that represents Baryogenesis without Lorentz violation (or 
Vb < vx)- As the Baryogenesis temperature is very high (Grand Unification Level), the Lorentz violation 
effect, when the Baryogenesis starts, it is very important; whenever vt > vx- The effects of this factor 
on the solution will be seen in a subsequent work. So far, the important result is that it is possible to 
find a trace of a possible Lorentz violation in the Baryogenesis. 

Remembering the bound found with the Threshold Energy for the boson decay, if Vb > vx- 

Px < -f= (73) 
V2aa 

we can find a limit to the temperature when these reactions start. For this, we are looking for the 
temperature to fuUfill that: 

mxc 
v2oa 

For this, we need the relation between average momentum and temperature. Using a Fermi statistic 
and Ex = vxPx, we obtain: 

So, the temperature at the beginning of the Baryogenesis is: 

30C(3)mxc2 



^4 ~ 0.3702 X 10" (76) 
mxC' 

Where we used da = 5 x 10~^^ [TT]. As the energies are in the Grand Unification level, it is required 
that ksTs or nixc^ ^ 10^^ [GeV]. If we impose this limit to mx, we obtain a temperature: 

kBTs > 0.3702 X lO^'^iGeV] (77) 

That matches a too much early era in the universe (Planck era). On the other hand, if we impose the 
limit to Tg: 
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mxc^ > 2.7012 X 10^[GeV] 



(78) 



So, in spite of having an extremely high mass (tox 3> mi,), these values are far below of the Grand 
Unification level (Desert). So, it is possible that the X Boson would be observed in the LHC because the 
maximum energies are ^/s = 14 [TcV] in proton-proton collisions [13] ■ 



CONCLUSION. 

As we are at a high energy level, in the Grand Unification scale, we have a greater possibility to find 
a LIV effect. In this work, we saw that this effect really exists, but it becomes more important if the 
baryon and boson maximum velocities are related by Vb > vx- Owing to the fact that we do not know 
the X boson properties and difhculty to estimate the a parameter, it is difhcult to know if we are really 
in this case. But if this is the case, we would have an important trace of a LIV. Additionally, it is possible 
that the modification to the differential equation of the X boson decay will give some information about 
the probability amplitudes. It could be interesting if a LIV could be interpreted as a parameter of the 
type of e in (l6|). If this is possible, the LIV could explain the preference for the matter over the anti-matter. 

We estimated a condition for the moment when the Baryogenesis begun, given by the LIV. This 
condition tells us that ksTs = 0.262 x lO^^mxc^. Then the majority of bosons start to decay. The 
condition give us a estimation of the mass or temperature if we have other method to know one of them. 
Anyway, it is a condition that only appears if we impose the LIV. 

Other point to analyze, that we leave proposed, is the possibility that the bosons arc part of the dark 
matter. The inferior limit found to the boson mass is very high yet, but it is very close to the energy limit 
obtained by the LHC, so the X boson could be produced. We must say that this mass limit is obtained 
only by the LIV too. 
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